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A detailed analysis of the small-amplitude solutions of a deformed discrete nonlin- 



" singular" points which split the infinite chain in a number of independent segments. 
■ We show that small-amplitude dark solitons in the vicinity of the singular points 

o: 

are described by the Toda-lattice equation while away from the singular points are 
described by the Korteweg-de Vries equation. Depending on the value of the defor- 
mation parameter and of the background level several kinds of solutions are possible. 
In particular we delimit the regions in the parameter space in which dark solitons 
are stable in contrast with regions in which bright pulses on nonzero background are 
possible. On the boundaries of these regions we find that shock waves and rapidly 

X. 

^ . spreading solutions may exist. 
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I. INTRODUCTION 



The general discrete nonlinear Schrodinger equation (GDNLS) 

iq n + (1 + r]\<ln\ 2 ){qn~i + q n +i ~ 2<?n) + 2(u n + \q n \ 2 )q n = 0. (1) 

was introduced in ref. |T[ as a generalization of the simple tight-binding linear Schrodinger 
model for the dynamics of quasiparticles in a molecular crystal. In this equation q n rep- 
resents the complex mode amplitude of the molecular vibration at site n, u n is the onsite 
frequency of the vibration while the nonlinear terms arise, in the adiabatic and small field 
approximation, as the result of the interaction of the quasiparticle with the lattice. From 
a mathematical point of view this system represents a norm preserving deformation of the 
diagonal discretization of the nonlinear Schrodinger equation (DNLS). The presence of the 
deformation parameter rj G 9ft in the system allows one to study, both at the classical and at 
the quantum level, the interplay between onsite-intersite interactions as well as integrability- 
nonintegrability and discrete-continuum properties 0-^,3- For rj = 1 (off-diagonal nonlin- 
earity) the GDNLS reduces to the integrable Ablowitz-Ladik (AL) model with exact soliton 
solutions while for 77 = (onsite nonlinearity) it gives the nonintegrable DNLS system. For 
intermediate values of rj it allows all possible splitting of the nonlinearity among three adja- 
cent sites (equal splitting is at 77 = 2/3). The properties of the localized mode solutions and 
of the plane waves (modulational instability) of eq.(l) were studied in a series of publications 
@||. Bloch oscillations of bright solitons in the external electric field (u n = ^n, 7 being a 
constant) at 77 = 1 were also considered in || and in the case < 77 < 1 in while the 
spatial properties and the existence of breather-like impurity modes were reported in refs. 
||, 0. All these studies dealt mainly with the attractive case (i.e. rj > 0) at uj n = 0, in 
which brigth solitary pulses are stable. 

The aim of the present paper is to investigate the properties of the small amplitude 
solution of the GDNLS equation in the repulsive case rj < 0. In the following we refer to this 
as the stable case since the plane wave solutions corresponding to the center of the Brillouin 
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zone (BZ) are modulationally stable. To this end we consider rj = —e with < e < 1 and 
fix u n = p 2 with p a constant which will be associated with the amplitude of q n either at 
the infinity or at the singular points (see below). 
The GDNLS is then written as 

iq n + (1 - e|g n | 2 )(g n _i + q n+1 - 2q n ) + 2(p 2 - \q n \ 2 )q n = 0. (2) 

The inverse scattering technique for (§) at e = 1 (stable AL) has been developed in HlU| , its 



dark soliton solution has been found in [I3J and Bloch oscillations of the dark soliton in a 
constant electric field were reported in [12| (note that we have changed the sign in front to 
the onsite nonlinearity in eq.^ just to have e = 1 as the AL limit). 

Though many aspects of the behavior of solutions of (|2]) at e = 1 are very similar to the 
behavior of their counterparts in the continuum limit (see e.g. |13|]) it has been shown in 
|14] , p!5[| that there are drastic differences between the discrete and continuum dynamics when 
the amplitude of a site excitation approaches one. In this case the dispersive term becomes 
zero and the given site decouples from its neighbors. It is a direct consequence of the explicit 
form of (0) that at e < 1 a site must have an amplitude p = e -1 / 2 to become decoupled. 
When this happens we call the corresponding point a singular point. Thus singular points, 
splitting the infinite chain in a number of independent segments, exist for the GDNLS for 
arbitrary nonzero values of the deformation parameter. In the present paper, by performing 
a multiple scale expansion of the GDNLS around singular points, we show the existence 
of dark solitons which (in the small-amplitude approximation) are described by the Toda- 
lattice equation. Away from singular points we find that dark solitons can exist in certain 
regions of parameter space (e, p) and in the small-amplitude limit are described by solitons of 
the Korteweg-de Vries (KdV) equation. In addition we find a region in the parameter space 
in which dark solitons are unstable and bright pulses on nonzero background are stable. On 
one of the boundaries of this region we find, quite surprisingly, that the system becomes 
effectively dispersionless and the formation of shock waves becomes possible. The second 
boundary corresponds to effectively linear regime where initially localized pulse spreads out 



quite rapidly. The analytical results, derived in terms of a multiple scale expansion, are 
found in good agreement with direct numerical simulations of the GDNLS even behind the 
limit of validity of the small amplitude approximation. 

The paper is organized as follows. In section II we study singular point solutions of the 
GDNLS by considering the case of a single site between two singular points. In section III we 
examine small-amplitude dark solitons in the neighborhood of the singular points in terms 
of the Toda chain. In section IV we consider a multiple scale expansion of the dynamics 
of chain excitations when the amplitude of the background is less than the amplitude of 
the singular points. Finally we briefly discuss the regime of zero effective nonlinearity and 
the possibility of shock wave solutions in the GDNLS. In the conclusions we summarize the 
main results of the paper. 



where k denotes the wave vector from the first Brillouin zone (BZ): k 6 [— 7T, 7r], uj{k) = 
4(1 — ep 2 ) sin 2 (£;/2) is the frequency of the background and $ is a constant phase. It possess 
the following integral of motion 



II. SINGULAR POINTS AND ONE-SITE DYNAMICS 



Let us consider equation (§) subject to the non-zero boundary conditions 



pe 
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As mentioned above, if at some site n the amplitude of oscillations is e 



/ 2 , and hence 



q n = 



— 1/2 iuJot— ikn+i'd 



(5) 



with 




(6) 



the evolution of g ra _i becomes disconnected from q n+ \. 
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So, for two singular points placed on sites /_ and /+ (being l± integers), there is an 
integral of motion 

i+ 

= E {e(q n q n -i + q n -iq n ) + 2(1 - e)|g n | 2 } (7) 

n=l-+l 

associated with the dynamics of the segment between l_ and /+, which naturally follows 
from (f|). 

The integral allows us to find an explicit solution in the case in which one point is 
placed between two singular points. For the sake of definiteness let us assume that l± = ±1, 
with 

q ±1 = e~ 1/2 exp(iuj t ± ztf) (8) 

being singular points. Representing the solution of the middle point as 

q = €~ x l 2 v exp(iuo t) (9) 

one arrives to the following expression for the integral of motion 

2_i,i = 2C cos ^ + i^i(l + C 2 + ^). (10) 

Here ( and £ are real and imaginary parts of v. v = ( + Just form this expression 
one can see an essential qualitative difference in dynamics of the AL model and GDNLS 
equation. In the first case ||15|| , i.e. at e = 1, there is a region of parameters (at |z/| > 1) 
where the lattice becomes unstable, this being displayed by an infinite growth of £ during 
finite time. In contrast, even a small value of the difference 1 — e will prevent this growth. 
The equations for (CO can be obtained from Hamilton's equations with respect to the 
nonstandard Poisson brackets 

t/if} _ [1 _ (< . + fl 1 {§£|_g|} 

with respect to the Hamiltonian 2_ 11 . Indeed, for 22 = I-i i we have 

C = {22,C} = -2^[l-(C 2 + a]£ (12) 
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i = {h, o = - (c 2 + em + o (is) 

where A = -r 1 - cos 

Taking into account the explicit form of one finds a convenient parametrization of 
the problem as 

C + A = i?cosx, £ = Rsm X , (14) 

where R is a positive constant playing the role of the radius of the orbit in the phase portrait 
and is related to the energy by 

R 2 = + A 2 - 1. (15) 



The dynamical equations (]T2|), ( |T3| ) are easily solved and give the following two different 
types of orbits. 

(a) If R < |1 — A| or R > |1 + A| (hereafter without restriction of generality we take 
4> £ [— tt/2, 7r/2]) the motion is periodic and is governed by the equations 

cos — 1 . . 

sin 



where 



9 = 4\A4 2 - LRcostf • t + 2arctan — A — ^ — - (18) 

\A + IR + A + ( J 

and 

A = 1 ~ R ] - A ' (A 2 >1). (19) 
2Ai? v 7 v 7 

Here £o> Co denote the initial values at t = 0. 

(b) If |1 — A | < i? < |1 + A | the system displays aperiodic motion described by the 
equations 



C = -R 



A cosh 6-1 



+ A 



(20) 



cosh 6 + A 



f = Ry/l - A 2 



sinh B 



(21) 



cosh 6 + A 



where 




(22) 



and the constant A is given by ( |T§D (note that now A 2 < 1). 



The existence of two regimes oscillating and aperiodic is evident from the phase portrait 
which consists of circles of the radius R which in the integrable limit, e = 1, degenerates 
into a straight line. In the aperiodic case a circular trajectory always crosses the unit circle 
(given by \u\ = 1), the crossing point being just the singular point. Thus, for \A\ < 1 an 
initial condition will move on the circle characterized by R and, as t — > oo, it will reach the 
singular point on the unit circle. For \A\ > 1 there will be no crossing and the point will 
continue to rotate on the respective circle. 

From the above analysis it follows that if initially \q n \ < e _1 (or \q n | > e" 1 ) this will be 
true for all times. This property is valid also for the integrable AL system. 

III. SMALL- AMPLITUDE DARK SOLITONS NEAR SINGULAR POINTS 

Rather complete description of the dynamics of (0) is also available for the small am- 
plitude pulses. Indeed starting with the case of the excitations slightly deviating from the 
singular points, we can consider backgrounds for the nonlinear excitations in the from of 
arbitrary (i.e. with arbitrary k values in the BZ) set of the singular points. For the sake of 
simplicity in the present section we restrict the analysis either to the center of BZ, k = 0, 
or to the edge of BZ, k = ir. 

The AL dark soliton (e = 1) is written in the form 



sin 2 $ 



n-i/2 



= P 1 



cosh 2 (Kn - tit) 



exp [— itan^tanh^n — fit)] 



(23) 



S 



with the parameters K, p linked by 



K p 

sinh — = . smfi. (24) 

2 v/T^P 7 



The phase •& parametrizes the family of the soliton solutions and is chosen in the interval 
[— 7r, 71"]. The small amplitude limit corresponds to the case •& <C 0. In order to expand ( |23|) 
around *& = we observe that the most natural representation of the solution of the GDNLS 
equation at e < 1 is 

q n = K n e~ 1/2 (l - 7 2 /ia„)e- iwx ' l+ ^ ' (25) 
where 7 1 is a small parameter, a n = a n (r) and \ n = x«( r ) are t wo real functions of 



the slow time r = 2 / yy2(e~ 1 — 1 — K)t, and luo is given by (|g). The parameter k — ±1 
introduced in ( Pq) will be used to choose two different background oscillations i.e. in-phase 
oscillations (k = 1) and out-of-phase oscillations (k = —1) corresponding to k = and 
k — 7T, while the parameter = ±1 is used to obtain either dark (p = 1) or bright (p = —1) 
pulses. 

Substituting (|25|) in (fj) and gathering all the terms of orders up to 7 s we arrive to the 
following Toda system |13|] 

da r . 



dr 

dc n 
dr 



(26) 

— a n+l — On; (27) 



where 



~n= \72le- 1 - l-K|(Xn-Xn+l), (28) 



and 



= Ksign(e 1 — 1 — ft). (29) 

(For the integrable case this system was also obtained in |TB[). We have therefore that small 
amplitude dark pulses near the singular points can be viewed as exact solitons of the Toda 



lattice moving on the oscillating background pexp(iuot). When k = 1 (k = 0), it follows 
from ( |29|) that \i — 1 i.e. the respective solution (|25| ) is always dark, i.e. its energy profile is 
always below the background level. On the other hand, at the edge of the BZ we have that 
k = —1 and the sign of \i depends on e. We find that \i — 1 for e < 1/2 and fi = — 1 for 
e > 1/2. This implies that at k = it small amplitude dark solitons exist for e > 1/2 while 
for e > 1/2 bright pulse should appear. To check these predictions we have numerically 
integrated the GDNLS system with a fifth-order Runge-Kutta adaptive stepsize algorithm 
(the numerical error was controlled by checking the conservation of the norm up to the fifth 
decimal digit). In Fig. 1 the time evolution of a small amplitude dark soliton on a line of 
1500 sites with an in-phase (k = 1) background is reported. Here the initial condition is 
of the form of a one-soliton solution of the Toda lattice which corresponds to a bisoliton 
solution of the GDNLS. We see that the initial pulse splits in two dark excitations which 
are stable over long time. This is in contrast with what happens for a bright profile on the 
same background and for the same parameter values as reported in Fig.l. This behavior 
is found to be true for all values of e in the range (0, 1]. The same analysis, but for an 
out-of-phase background (k = —1), shows that bright pulses are stable and dark ones are 
unstable if e < 1/2, while dark pulses are stable and bright ones are unstable if e > 1/2. 
In Figs. 2,3, we have reported the evolution of a bright and dark GDNLS bisoliton for 
e values respectively at e = 0.25 and e = 0.75. We see that the pulse are quite stable in 
agreement with our predictions. From the above analysis it is also clear that at the " critical" 
value e = 1/2 the solution changes stability and the dynamics manifests a strong dispersive 
character (weakly nonlinear or effectively linear regime). In Figs. 4,5 the time evolution of 
a bisoliton, respectively dark and bright, at the "critical" value e = 1/2 is reported. We 
see that in both cases the profiles decay in background radiation in agreement with our 
prediction. 
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IV. SMALL- AMPLITUDE EXCITATIONS FAR FROM THE SINGULAR POINTS 

AND SHOCK WAVES OF THE GDLNS 



In order to find dynamical equations for small-amplitude pulses far from the singular 
points we apply the same analysis as in the previous section but now we expand around 

q n =[p + a n (t)}e- l ^ t \ (30) 

where 

a n (t) = 7 V°)(T,X;r) + 7 V 1 )(T,X;r) + ... (31) 

$ n (t) = 7$ (0) (^ X; r) + 7 3 $«(T, X; r) + ... (32) 

In the above expansion X = ^n, T = jt denote fast space and time variables while 
t = 7 3 t represents a slow time. By analyzing the equations of all orders from 7 up to 7 s we 
find that at the leading orders (7 2 , 7 s ) the GDNLS system reduces to 

«9$(o) 



dT 



4pV°J (33) 



-w = (1 - e ^w (34) 



From system (|33D , (0) it readily follows that the excitation at leading order moves with 
the velocity 



V e = -2p^fl-ep 2 (35) 

[the given sign of the velocity is chosen only for the sake of convenience] which depends on 
the background level. 

It is of interest to note that the velocity V = Q/K of the AL dark soliton in the small 
amplitude limit, takes the form 

v = v °+ of P 2 (3~V), (36) 
11 



where Vq is given by (|35|) at e = 1 and coincides with the velocity of linear waves against 
the background p in the long wavelength limit (i.e. in the center of BZ). It also follows from 
(|36|) that V coincides with the velocity of the linear waves when p 2 = 3/4. 
In the next two orders in 7 we obtain the equations 

-fffr - Va (1) = --^r - (1 - ep 2 ) w + 6pV°) 2 , (37) 



<9T v ^ ; dX 



2 



^ + l^_ (2 _ 5£p2) (^ o( .,^. (38) 



dr 12 <9T<9X 2 v r ' V V e J dX 
It is remarkable that the compatibility condition of these two equations directly leads to 
the following KdV equation 



- 4 Pv /l - ep 2 ^ - -(1 - ep 2 )[3 - (3e + l)p 2 ]^F + 8 -° 2 ( 3 - 4e P> ( ° ] = °> ( 39 ) 

where Z denotes the running variable Z = X — V t T. We remark that a similar result was 
also obtained in refs. ||. 

Equation ( p9|) allows soliton solutions corresponding both to positive and negative 
depending on the sign of the prefactors in the dispersive and nonlinear terms. The results of 
the respective analysis are illustrated in Fig.6. In this figure the region between the curves 1 
and 2 marked with B corresponds to parameter values for which stable propagation of bright 
solitons against the background is possible. In this region dark solitons will not be stable 
and they will decay in background radiation. The line 3 itself determines the amplitude of 
the singular points: e = 1/p 2 . Close to this line the dynamics of small amplitude excitations 
follows Toda equations as described in section 2. The areas marked with D in Fig. 6 
correspond to parameter regions in which dark solitons are stable while bright ones are 
unstable. On the curve 2 of Fig. 6 [e = 3/(4p 2 )] the nonlinear term becomes zero as it 
is evident from (|39|) . This implies that any localized pulse will spread out in background 
radiation (no soliton-like excitations). On the other hand on curve 1 of Fig. 6, which 
represents the other dark-bright interface, the system has different dynamical properties 
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since it becomes effectively dispersionless as it is seen from fl3~9"|). Since the nonlinearity on 
this curve is not zero we expect this case to give shock solutions (see below). 

To check these results we have numerically computed the time evolution of small am- 
plitude bright and dark excitations for parameter values taken in different regions of Fig. 
6. In Fig. 7a the evolution of a dark soliton at e = 0.4, p = 0.8 is reported. In Fig. 7b 
we show the same evolution but for an initially bright pulse for the same parameter values 
as in Fig. 7a. From these figures it is clear that parameter values on the left of curve 1 
correspond to stable dark pulse propagation and to decaying bright excitations. A similar 
behavior is found in the region between curve 2 and curve 3. This is seen from Figs, 8a,b in 
which the evolution of, respectively, a dark and a bright excitation is reported at parameter 
values e = 0.4, p = 1.5. In the region between curve 1 and curve 2 the behavior is just the 
opposite i.e. dark solitons decay in background radiation while bright pulses can propagate 
as solitary waves. This is shown in Figs. 9 a,b for respectively a bright and a dark initial 
profile at parameter values e = 0.16, p = 1.7. It is of interest to remark that while for the 
pulses of Fig. 7b, the radiation is on the front in Figs. 8b, 9b the radiation is just on the 
back. This corresponds to the fact that wave packets in the region on the left of the curve 
1 always have group velocity greater than the velocity V e of the background radiation while 
the opposite is true for the complementary regions. 

From (|39"D it is evident that for e = 1/p 2 — 1/3 the dispersion disappears and the resulting 
equation possesses shock wave solutions. This surprising property of the GDNLS system 
has been numerically checked in Fig. 10. In this figure the time evolution of an initial 
bright profile is reported for parameter values taken on the curve 1 of Fig. 6 at e = 0.361(1), 
p = 1.2. We see that the smooth initial profile is distorted with the top moving at higher 
velocity of the bottom part of the profile. This give rise to a forward bending of the profile 
with breaking of the wave and rapid oscillations on the wave front. A similar behavior is 
observed also for dark pulses as it will be reported in more extended form elsewhere J20| . 



We should mention that the possibility of shock waves in other chains was also reported in 
ref. [gl . 
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V. CONCLUSION 



We have analyzed the properties of the small-amplitude solutions of the GDNLS. Like the 
well-known AL model the GDNLS equation has singular points for all e different from zero 
which provide effective decoupling of the chain. In this sense one can speak about dominant 
role of the nonlinearity of the "AL-type" in the behavior of the system. Meantime, the 
properties of the singular points are essentially different in the cases e = 1 and e < 1. 

The small amplitude excitations against constant amplitude backgrounds are governed 
by the Toda-lattice equation if they are in the vicinity of the singular points or by the 
KdV equation if the background amplitude is less than 1/y/e. The type of the solution 
essentially depends not only on the deformation parameter but also on the chosen point in 
the BZ. Existence of the bright pulses against nonzero background above singular points is 
a nontrivial qualitative demonstration of the mentioned differences. We have delimited the 
regions in the parameter space in which dark solitons are stable in contrast with regions in 
which bright pulses on nonzero background are. On the boundaries of these regions we have 
found, quite surprisingly, that shock waves may exist. 
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FIGURES 

FIG. 1. Time evolution of an initial dark pulse (p = 1) on an in-phase background (k = 1) at 
e = 0.25. The total integration time is 300. The initial condition corresponds to a soliton of the 
Toda lattice. 

FIG. 2. Time evolution of an initial bright (// = —1) pulse on an out-of-phase background 
(k = —1) at e = 0.25. The total integration time is 300. The initial condition corresponds to a 
soliton of the Toda lattice. 

FIG. 3. Same as in Fig.2 but for a dark initial condition at e = 0.75. 

FIG. 4. Time evolution of an initial dark pulse (p = 1) on an out-of-phase background 
(k = —1) at the critical value e = 0.5. The total integration time is 1200. 

FIG. 5. Same as in Fig.4 but for a bright pulse (p = —1). 

FIG. 6. Parameter space (e, p). The window where bright pulses against a nonzero background 
exist is bounded by the curves (1) e = l/p 2 - 1/3 and (2) e = 3/(4p 2 ). The curve (3), e = l/p 2 , 
corresponds to the singular points and in its vicinity the small-amplitude excitations of the system 
are governed by the Toda-lattice equation. 

FIG. 7. (a) Time evolution of an initial dark pulse corresponding to a KdV soliton for 
parameter values e = 0.4, p = 0.8. (b) Same as in (a) but for initial bright pulse. The total 
integration time is T = 1200. 

FIG. 8. (a) Time evolution of an initial dark pulse for parameter values e = 0.4, p = 1.5. (b) 
Same as in (a) but for an initial bright pulse. The total integration time is T = 1200. 

FIG. 9. (a) Time evolution of an initial bright pulse for parameter values e = 0.16, p = 1.7. (b) 
Same as in (a) but for an initial dark pulse. The total integration time is T = 1200. 
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FIG. 10. Shock wave formation from a smooth bright pulse at parameter values 
e = 0.361(1), p = 1.2 on curve 1 of Fig. 6. The total integration time is T = 1200. 
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